Controlling quantum systems by embedded dynamical decoupling schemes 
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A dynamical decoupling method is presented which is based on embedding a deterministic de- 
coupling scheme into a stochastic one. This way it is possible to combine the advantages of both 
methods and to increase the suppression of undesired perturbations of quantum systems signifi- 
cantly even for long interaction times. As a first application the stabilisation of a quantum memory 
is discussed which is perturbed by one-and two-qubit interactions. 
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The desire to control quantum systems and to push 
characteristic quantum phenomena more and more into 
the macroscopic domain necessitates the development of 
suitable methods which are capable of cancelling unde- 
sired perturbations as much as possible. Such methods 
are not only of significant interest for the further devel- 
opment of quantum information science but are also im- 
portant in more traditional areas of quantum physics, 
such as high-resolution spectroscopy. So far, two differ- 
ent strategies have been pursued to suppress undesired 
perturbations in quantum systems. First of all, one may 
cancel such perturbations by appropriate sequences of 
control measurements and unitary recovery operations. 
This approach is at the heart of quantum error correc- 
tion m. it requires the use of redundant quantum states 
and, in principle, it is capable of correcting errors per- 
fectly. However, with nowadays experimental techniques 
it is difficult to realize typical syndrome measurements 
and recovery operations required. 

The strategy pursued in dynamical decoupling meth- 
ods [ll!llllla0&HE[l3 is different. Thereby, one 
uses controlled unitary dynamics to suppress undesired 
perturbations as much as possible. As these methods are 
not based on redundancy and quantum measurements 
they tend to be implementable in realistic scenarios more 
easily. Depending on the type of dynamics applied 
one distinguishes deterministic and stochastic decoupling 
methods. Deterministic methods HHSHISIlllla 
are already well developed. Numerous deterministic se- 
quences of unitary operations are already known which 
are able to suppress few-qubit perturbations. Stochas- 
tic decoupling methods [lfjt lllj have been introduced 
only recently. In particular, it was demonstrated that for 
sufficiently long interaction times stochastic decoupling 
methods offer advantages over deterministic ones. Typi- 
cally, they lead to a linear-in-time fidelity decay whereas 
the fidelity decay of deterministic methods is quadratic 
in time. Nevertheless, for sufficiently short interaction 
times the error suppression capabilities of deterministic 
methods are typically stronger than of stochastic ones. 

In this letter it is demonstrated that the advantages of 
both methods can be combined by embedding an appro- 
priate deterministic decoupling scheme into a stochastic 
one. As a result, one achieves the significant suppression 
of perturbations which is characteristic for deterministic 



decoupling schemes and, in addition, one also produces 
a rather slow linear-in-time fidelity decay which is typi- 
cal for stochastic decoupling methods. Motivated by the 
potential impact these error suppression methods might 
have on the future development of quantum information 
science in the following the basic ideas of these embed- 
ded decoupling methods are demonstrated for a quantum 
memory which is perturbed by inter-qubit couplings. 

Let us start by summarising briefly the main ideas 
underlying deterministic dynamical decoupling methods 
HH H IE H H H in more detail. For this purpose 
we consider a quantum state \ip) of a n g -qubit quantum 
memory which is perturbed by an unknown inter-qubit 
coupling Hamiltonian, say Hq. The resulting coherent 
errors may be suppressed by modifying the n 9 -qubit dy- 
namics with the help of an appropriate periodic deter- 
ministic control sequence of unitary operations. Within 
the framework of the bang-bang methods[2(, for example, 
such a periodic cycle of duration T c consists of N unitary 
operations. They are of the form djdj_ 1 with properly 

chosen unitary operations dj which are applied (approx- 
imately) instantaneously at times to,...,ijV-i (compare 
with FigUJi). Thus, after one cycle the resulting dynam- 
ics of the quantum memory are described by the unitary 
time evolution operator |10| 
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with to = 0, = T c , and with the 'interaction-picture' 
Hamiltonians Hj = d^Hodj ■ In order to suppress the per- 
turbing influence of Hq as much as possible this control 

cycle should fulfil the requirement X^o 1 -^j'G^j+i ~ 
0. This condition guarantees that at multiples of the 
cycle time the perturbation is cancelled by the control 
sequence at least to lowest order perturbation theory. 
Nevertheless, as this cancellation is not perfect, repeated 
applications of the deterministic control cycle lead to a 
coherent accumulation of the residual perturbation char- 
acterised by the Hamiltonian H with 
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(a) didj dodjv-i 
-^H 1 

time tN+1 t N = T c 
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(c) didj dofifjdjy-i 
-<H 1 

time tjv+1 tjv = Tc 



FIG. 1: Schematic representation of the deterministic bang- 
bang method (a), the PAREC- method (b), and an associated 
embedded decoupling scheme: The time evolution between 
subsequent instantaneously applied unitary operations is gov- 
erned by the perturbing Hamiltonian Ho. In the embedded 
decoupling scheme (c) the deterministic decoupling cycle of 
(a) is embedded between any two subsequent random unitary 
operations of the PAREC- method (b). 



As a result any linear superposition of eigenstates of H is 
dephased. Thus, the short-time evolution of the fidelity 
f(T) of the quantum state 1-0} is given by [T^j 



f(T) 
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with the relevant energy uncertainty of the residual inter- 
action (AH) 2 = {il>\{H- (-)p\H\ip)) 2 \ip). An upper bound 
of the residual interaction is given by the norm of its 

Hamiltonian \\H\\. Using Eq.J5J) and basic properties of 
the logarithmic function we may obtain the inequalities 
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Thus, for small cycle times T c with 
lower bound of the fidelity is given by 
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In general, deterministic decoupling schemes suffer from 
two major drawbacks Firstly, some perturbations 

may require rather long control cycles so that a signifi- 
cant error suppression cannot be achieved. Secondly, it is 
unclear how to suppress perturbations which change on 
time scales short in comparison with cycle times required. 

In order to overcome these drawbacks stochastic de- 
coupling methods have been proposed recently 0, ^| . 
In the Pauli-random-error-correction-(PAREC)-method 
[Tlj . for example, the periodically repeated control cycles 
of deterministic methods are replaced by a single random 



sequence of (approximately) instantaneous uncorrelated 
unitary operations. Thus, the deterministic unitary op- 
erations dj of the bang-bang methods are replaced by 
unitary transformations fj each of which is a product 
of uncorrelated random one-qubit Pauli-operators, i. e. 
fj 6 {i, X ,Y , Z}® n " , with the Pauli spin operators X, 
Y, and Z. Within a time interval of duration T the re- 
sulting dynamics are described by a unitary time evolu- 
tion operator of the form of Eq. with the replacement 
T c — > T and with the 'interaction-picture' Hamiltonians 

Hj — > fjHofj. In order to suppress the perturbing influ- 
ence of Hq the random sequence of unitary Pauli oper- 
ations fj should be statistically independent. Provided 
these uncorrelated random Pauli-operations are applied 
after equal time intervals of duration At a lower bound 
of the mean fidelity is given by 0, 0, 0] 
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if the interaction time T = nAt is so small that | 
I — E/(T) |« 1. (E denotes the statistical averaging 
over the uncorrelated random unitary operations fj.) 
Consistent with Fermi's Golden rule the decay rate is 
given by V = (2-K/h)[At/ (2-kK)]AH 2 and (AH Q ) 2 = 
(4>\(H ~ (tp\H \ip)) 2 \ip) is the relevant energy uncertainty 
of the original perturbation. The lower bound of © 
is based on the inequality (AH ) 2 < ||H IP- 
According to relation © stochastic decoupling meth- 
ods, such as the PAREC-method, offer the advantage 
that residual errors which are not corrected by the ran- 
domly applied unitary operations do not add coherently 
so that the fidelity decays linearly with interaction time 
T and not quadratically as in the case of deterministic 
methods. As a consequence, for sufficiently long interac- 
tion times stochastic decoupling methods offer a higher 
degree of error suppression than deterministic ones. But, 
for sufficiently short interaction times deterministic de- 
coupling methods are still expected to be more accu- 
rate because the lower bounds of their fidelities scale 
with ||-ffo|| 4 (compare with (J5J) in contrast to the ||7?o|| 2 - 
scaling predicted by relation JJJJ. 

By embedding a deterministic decoupling scheme into 
a stochastic one errors are suppressed to a large extent 
not only for short but also for long interaction times with- 
out requiring any major additional effort. The basic idea 
of such an embedding scheme based on the recently pro- 
posed PAREC-method is depicted in Fig^. In order to 
stabilise a quantum state over a time interval of mag- 
nitude T we first of all take advantage of the PAREC- 
method by applying (approximately) instantaneous ran- 
dom Pauli operations after each time interval of duration 
At. In addition, in order to get rid of large parts of the 
undesired perturbation Hq within each of these time in- 
tervals we embed a deterministic decoupling scheme with 
cycle time T c — At between any two of these random 
Pauli operations. As a result, within each time interval 
of duration At the time evolution is now governed by the 
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residual interaction H . Furthermore, large parts of it are 
suppressed by the randomisation of the PAREC-method. 
Thus, this embedding procedure decreases the fidelity de- 
cays of Eqs.© and © significantly. In particular, using 

the replacements Hq — > H, At — > T c in the lower bound 
of Eq.©, and invoking relation Q for \\H \\T c /h < 1 
one obtains the lower bound 
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for the mean fidelity of the resulting quantum state. 
It indicates already a significantly slower fidelity decay 
than the corresponding decays of the deterministic and 
stochastic procedures of relations JSJ) and JBJ. In partic- 
ular, this lower bound of the fidelity decays linearly with 
T and, in addition, it scales with ||i?o|| 4 - 

Let us now investigate this fidelity decay quantitatively 
in more detail. For this purpose we consider a quan- 
tum memory consisting of n q — 9 distinguishable qubits 
which are arranged in an equidistant two-dimensional ar- 
ray (compare with FigEJ - This quantum memory is as- 
sumed to be perturbed by Ising- and Heisenberg-type 
qubit interactions of the form Hq/H — Ylk=o ^h^k + 
Y,k<i=o Jki(X k Xi+Y k Yi + Z k Zi) with the time indepen- 
dent perturbation parameters S k and J k i- Initially, the 
quantum memory is assumed to be prepared in the coher- 



ent state 
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with the dimension of the Hilbert space D = 2 nq and 

with |m = = l*i> ® l*a) ® • • ' \in g ) in £ {0, 1}) 

denoting the D orthonormal basis states of the compu- 
tational basis. 

The resulting time evolution of the fidelity is depicted 
by the full lower curve of FigEl F° r sufficiently short 
interaction times T the fidelity decay is given by 

f(T) = l-(TAH /h) 2 +G((TAH /H) 3 ). (8) 

The resulting quadratic-in-time decay reflects the de- 
phasing of the initially prepared quantum state \ip). At 
larger interaction times the dynamics are governed by 
recurrence- and revival-phenomena which characterise 
the intricate quantum interferences involved |14|. The 
extent to which the PAREC-method is capable of stabil- 
ising the initially prepared quantum state is exemplified 
by the full upper curve of Fig[3 The fidelity was obtained 
by averaging over 200 numerical runs. As expected, the 
resulting mean fidelity decay is linear in time and is well 
described by the decay rate T of relation JBJl (dashed up- 
per curve of Fig|2Jl. 

As the perturbing interaction Hq involves one-and two- 
qubit operations only an efficient deterministic decou- 
pling cycle can be realized by 32 deterministic unitary 
operations dj applied at equidistant time steps r within 
the framework of a bang-bang decoupling scheme H, || . 
These operations are products of Pauli-spin operators of 
the form d = cr^ ® • • • <& cr/ n with l k 6 {0, 1, 2, 3} and 
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FIG. 2: Time evolution of the fidelity f(T): without error cor- 
rection (full lower curve), resulting from the PAREC-method 
(full upper curve). The dashed curves indicate the approxi- 
mations of @ and JSJ. The interaction time T is plotted in 
units of the time interval r between subsequent randomly ap- 
plied uncorrelated Pauli-operations. The rectangular array of 
the nine qubits constituting the quantum memory is depicted 
in the inset. The lines indicate the non-zero values of the cou- 
pling parameters Jkii~. They are chosen randomly from the 
interval [— V3 x 10~ 3 ,\/3 x 1(F 3 ] (static imperfection). 
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FIG. 3: Fidelity decays resulting from the deterministic bang- 
bang method (full middle curve), the PAREC-method (full 
lower curve), and from the associated embedded decoupling 
scheme (full upper curve): The corresponding approximate fi- 
delities are indicated by dashed lines. The interaction time T 
is plotted in units of the time interval r between subsequently 
applied unitary operations. The inset shows the correspond- 
ing double- logarithmic plot. 



with So = 1,5-1 = X,& 2 = Y,a 3 = Z. The 32 rele- 
vant values (Zi, l nq ) are given by the orthogonal array 
OA(32,9,4,2) of Ref.[l||. The resulting time evolution 
of the fidelity is depicted by the full middle curve of FigOl 
It is apparent that this deterministic decoupling method 
leads to a significant slowing down of the fidelity decay 
in comparison with the uncorrected time evolution (full 
lower curve of Fig[3J). Furthermore, for the interaction 
times considered this deterministic decoupling procedure 
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also achieves a higher degree of stabilisation than the 
PAREC-method (lowest full curve of Fig El . This may 
be traced back to the rather small residual interaction. 
Furthermore, this fidelity decay is well approximated by 
relation Nevertheless, in view of its quadratic- in- 

time decay for sufficiently long interaction times T this 
deterministic decoupling scheme will eventually become 
less effective than the linear-in-time decaying PAREC- 
method. 

The fidelity decay originating from the embedded de- 
coupling scheme in which this deterministic scheme is 
embedded into the PAREC-method is exemplified by the 
full uppermost curve of Fig|3| Again the fidelity was ob- 
tained by averaging over 200 numerical runs. Over the 
whole range of interaction times the resulting fidelity is 
significantly higher than the corresponding values of both 
the deterministic and the stochastic decoupling method. 
Thereby, the lower bound of inequality J7J) strongly un- 
derestimates the actually achievable degree of stabilisa- 
tion. A much better approximation of the relevant fi- 
delity decay is obtained by evaluating the rate T of Eq. © 

on the basis of the residual interaction H. The slow 
linear-in-time exponential decay of the embedded scheme 



and the more rapid quadratic-in-time exponential decay 
of the deterministic scheme are also apparent from the 
inset of FigEl 

In summary, a generalised decoupling method has been 
presented which is based on the embedding of a deter- 
ministic scheme into a stochastic one. Thus it is possible 
to combine the strong error suppression properties of de- 
terministic methods at sufficiently short interaction times 
with the particularly advantageous properties of stochas- 
tic decoupling methods at longer interaction times. The 
example presented demonstrates the high degree of er- 
ror suppression achievable by such embedded decoupling 
schemes. Though this example has focused on the sta- 
bilisation of a quantum memory's quantum state, it is 
expected that embedded dynamical decoupling schemes 
also offer new and promising perspectives for stabilising 
the dynamical evolution of quantum systems or of quan- 
tum algorithms against undesired perturbations at least 
in cases in which the relevant deterministic method is 
able to achieve this goal. 
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